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ABSTRACT 


The object of this report is to investigate the response of a shaft 
or beam to shock applied transversely at the end supports. As a secondary 
objective, a simplified method of designing a shaft for shock resistance 
is suggested. In particular, the results of this work demonstrate the 
prediction of the transient history, with respect to time and position, 
of the beam deflection, shear force, and bending moment and show a calcit- 
lation of the times, positions, and maximum amplitudes of these quantities. 


To accomplish the objective, velocity shock was assumed and the 
differential equation for flexural vibration of Timoshenko was solved. 
the solution was then placed in a form such that values depend only on 
the ratio of the beam length to radius of gyration of a cross section and 
the ratio of shear modulus, modified by a factor accounting for section 
Shape, to the modulus of elasticity. The maxima for deflection, shear and 
bending moment were then tabulated for various length to radius of gyra- 
tion ratios. 


As a model, only a solid, circular, steel shaft of six inch radius 
was used, but the method of solution could be extended to other config- 
urations. 


As a result of the investigation, it can be shown that, comparing all 
Simply supported uniform beams of constant ratios of Length to radius of 
gyration of a cross section and of constant ratio of shear modulus to 
Young's modulus, modified by section shape, the time histories of deflec- 
tion, shear, and bending moment will differ only in magnitudes, i.e., the 
histories of any such beam can be scaled from that of another. 


It was found that there is a range of length to radius of gyration 
ratios where rotary inertia of elemental cross sections and shear effects 
Deecenewclenlt+cant. Burther, there are a two fold infinity of natural 
frequencies, each of which corresponds to an odd numbered deflection mode. 





The higher series of frequencies is associated with axial compression 
and tension due to bending, and the lower series is associated with the 
shearing effect. It was found that the higher frequency series, for all 
modes, goes to infinity as the length to radius of gyration ratio goes 
to infinity and that, as a result, the equations derived reduce to the 
classical flexural vibration equations. 


It is believed that the results obtained in this work could provide 
a designer concerned about the effects of shock with a useful tool in 
anticipating design performance due to shock. 


There was an indication that the energy contained in a beam subjected 
to shock cannot be adequately described with considerations of overall 
deflection. However, this point should be investigated further. 


It is recommended that the equations and curves obtained be checked 
by experiment. It was noted that the equation for shear converged rather 
Slowly; so, even though the curves obtained from computer calculations 
included terms through the forty-ninth mode, this may not provide enough 
accuracy. It is further recommended that other types of supports, 
especially fixed ends and cantilever, be investigated. 


Thesis Supervisor: F. M. Lewis 


Ta virers Professor of Marine 
Engineering, Emeritus 
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I. INTRODUCTION 


Today, through use of such devices as armor plating, subdivision, 
and torpedo bulkhead systems, naval ships are well able to remain afloat 
after being struck by sizeable amounts of explosives. However, a ship is 
useless if the shock associated with an explosion causes vital machinery 
to become inoperable. Furthermore, with modern nuclear weapons, depth 
charges, and mines available, physical contact with the ship is not 
always possible; and shock effects are then the primary means of incapa- 
citating ships. This is particularly true for submarines. Indeed in any 
industrial plant or ship where shock or jolting is anticipated, a design 
must be able to tolerate these effects in order to perform its mission. 

Of particular susceptability to the effects of shock are shafts and 
beams when their bearings or supports accelerate perpendicular to the 
shaft axis. In actual situations, such items as piping, condenser shells, 
rotating shafts, and in some cases plates could fall into the general 
category of beams. Of interest to the engineer are such problems as the 
amount of deflection that must be anticipated in order to have sufficient 
clearance around the shaft, and the magnitude of the stresses involved 
in order to assure that no plastic deformation occurs. Thus, this inves- 
tigation attempts to predict these values. 

There have been many works on the dynamic flexural behavior of 
beams. Among the primary of these is that of Timoshenko (ale However , 
these works neglect the effects of shear and rotary inertia which becomes 


Significant when the length to radius of gyration ratio becomes small. 


Eee 





To the author's knowledge, there is no simple design method established 
to account for shock effects. 

The investigation was limited to a solid, circular beam, simply supported. 
The effects of spinning, such as in rotating machinery, on shaft response was 
considered negligible; and no calculations account for this. 

The excitation was considered to be a "velocity shock", i.e. the beam 
Supports experience an instantaneous constant velocity at the incidence of 
the shock. This is a commonly used way of specifying shock excitation. Use 
of velocity shock allows calculations to be simplified; for when 2 body, 
initially at rest, is excited by a velocity shock, its response relative to 
the point of excitation is precisely the same as if the body were initially 
in motion with a velocity equal to the excitation velocity, and suddenly 
stopped ae Thus calculations can be made as though the beam supports are 
stationary and an initial constant velocity exists along the length of the 
beam. 

If one looks at a beam subjected to shock in actual practice, the response 
would depend upon the direction from which the excitation eminated. Thus, if 
one Support were excited before the other, the response would be different 
than if both ends were excited at the same time. In an attempt to simplify 
difficult calculations, this report is limited to the case of symmetrical 
excitation. : 

Linearity has been assumed in utilizing the equations involved in this 
report. Thus, the results are limited to beams whose stresses always remain 
in the elastic region. It was felt that this would not place a serious 
handicap on a designer if it is assumed that a beam has failed when the stresses 
resulting from calculations demonstrated herein exceed the yield limit of the 
material involved. 


= < 





To sum up, then, the model used for the investigation is a solid circular 
beam, Simply supported, with ends fixed; it is subjected to a velocity shock 
along its entire length. The response of the model is limited to the elastic 
limits of the beam material, which was assumed to be steel. The object of 
the work is to describe the response of the beam model with respect to deflec- 
tion, shear, and bending moment; then to show to what limits the results apply 
to other beams than the model chosen; and lastly to propose a simplified method 


of calculating the deflection and stresses for other beams. 





II. PROCEDURE 


Definition of Problem. Whereas the physical situation specifies that, 
at time, t =o, the ends attain an instantaneous velocity, v, the problem 
shall be considered as one where, at t =o, the entire length of the beam 
attains a velocity, -v. The final equation for deflection, then, will 
be the sum of the values derived and vt. 

Figure I indicates the model investigated and the sign conventions. 
Motion is restricted to the x,y plane. The problem is to find the motion 


of the beam indicated in Figure I. 


Derivation of the Differential Equation of Motion. Assume that the 
deflection is due to two separate and separable effects, shear and bending. 


Then total deflection is 


a (1) 
where 

ves deprecviensaue co shear only 

vy, = deflection due to bending only. 


Look now at Figure I. The elemental mass is considered rigid. 
Consider that, for all shear deflections, there is no rotation of the 


element. Rotation is associated with bending deflection only. Thus 


3 
Tv OM ae 
>. Moments = Qdx = So + J— ax=0 , 
: ot ox 


where the last term represents rotary inertia of the elemental mass. Thus 


it can be seen that, for bending moments, the element always remains 


Die 





Figure I. 


Beam Model and Sign Convention 
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Oy: 
perpendicular to the slope of the deflection due to bending, eee, and 





Ox 
ay 
5 Simply represents the angular acceleration of the element. 
ot ox 
From this it is clear that 
3 
OM oN, (2) 
Ol Ey ee 
rte sox 
Equating forces in the y direction, 
2 
\ forces = - oe Gx em 2 = 0 
Ox 2 ? 
ot 
or 
2 
) ) 
<2 - -m =f. (3) 
ot 


Recalling that the maximum shear stress, 
t 
© 5 (max ) cate: 


where } represents the slope of deflection due to shear, then 


oy 
an 


© 5 (max) ~ Ox 


Although shear stress is not constant over an entire cross section, it can 


be shown that the average shear stress, Se is 


oy oy 
1 t S S 
25 Soe oS ae 


Ge 






1 
where k is a factor depending upon section shape [2] ° 


Thus 


oy 


Q=S_A=- GA = , (4) 


Continuing, it is well known that 


32 
Vy 
il = NE — a (3) 
Ox 





Combining equations (1), (2), (3), (4), and (5), one obtains 


ae ss ay Se 
a, - (ES ) ay taz + FF ae (6) 
ox ax ot at ot 


Equation (6) is the equation to be solved. Its derivation is 
presented in reference [| » but has been given here in slightly different 


form because use will be made of the equations involved in the derivation. 


pOluUti@m Of Equetiongof Motions Assume that the variables are separable. 


Then 


vp (x,t) = (t ) X (x) 2 


The equation in x is then of the form 


4 2 
d x qd xX 


where A, B, and C are functions of time. The solution of equations of this 
type is of the form 


x) = C, sink x+C, coskx+C, sinhkx+C, coshkx. (7) 


-[- 


Since the ends are pinned, deflection and bending moments are always 


zero at the ends. Thus the boundary conditions are 


x(0) = X.(0) + X(0) =0 

X(L) = X,(L) + X(L) =0 (8) 
ax, ax, 

oo = ao = 

dx Ix oe 


Putting equations (8) into (7) it can be shown that X (x) is of the 


Perm, 
X (x) = a C. Sin kK, cay (9) 


where kK = ce 3 n= PS eene ° 


see Appendix A for proof of this. 


Assuming no damping, a solution for T (t) can be of the form, 


ic) = Ge cos ® t + C. Sinw t , 
n n n n 


in e 
or 
fae) 
y(x,t) = S Can cos W t+ C, Sin w, + | sink x. (10) 
eal Si See 


Putting (10) into (6) it is easily seen that the frequency equation results, 


i.e., 
Im 4 EIm 2 2 h 
=z O, E + (== + x) © | wo + EI = 0 (11) 


fa 





Equation (11) shows that there are two frequencies associated with 
t= Kelle} kw? which fact will be amplified in a later section. Thus, the 
assumed ue is incomplete, since it does not account for both frequencies. 


Allowing for both frequencies, let 


© 
y(x,t) = = ke sin we +f Con cos wa + aon sin w,° 
N=1,3,5,-¢- 
+ Cy cos wnt | sink x. (12) 


From equation (12), it is seen that four initial conditions must be 
utilized to evaluate the constants. Thus, initially, the deflection of 
the beam is zero. From the definition of the problem, the initial 
velocity along the entire beam is the velocity of shock, -v. Lastly, 
due to rotary inertia, the slope and angular velocity of the deflection 
caused by bending cannot instantaneously change, and are thus zero [3] : 


Putting this in mathematical terms, 








Jy iae0)) = 0 (as3)) 
oy _ 
st -_ = -v (14) 
= = 0 (15) 
x 
: t=o 
a 
fo) Vp 
xot Pe. =a (16) 





Now, put equation (12) in (3) and integrate. Then 


2 | 
dQ 7 oy 7 : m 2 
| ae = =m | 2 dx = = a ls or sin seal 
N=1,3,5,-6- 


&) 
HI 


2 2 mates One 
+ a Con cos wt + Ws on Sin wt + Ws Cin c08 Wp, [cos k x 


+ Const. (17) 


But from Appendix A, it can be seen that, at x = Z Q = 0, thus const. = 0. 


Apply equation (4) to (1%), so that 


2 2 2 
iw) ’ a a r 
+ An Con cos Wat + WD, en sin W,° + Wo, Oi cos vont | sin kx 


(18 ) 


Again, boundary conditions preclude constants of integration. Thus; from 


equations (1) (12), and (18), 








0 De | 
- fo Cc. sin w,_t + C,. cos w__t 
Vp = GA 2 in Yin en vita 
nN=1255, 5s n | 
By 
+ aaa ~2n 
aak © Con sinw,,t + C, cos wot Sin kx. (19) 


Now the initial conditions can be applied. Applying equation (13) and 


looking at the nth tern, 


ae 





Ye (x,O0) = (Co. + Cy? sin kx = ©. 


a Con = eth : 


Applying equation (14) , 








oy. 
n 
cles | a Comore eceem 2 ~ -Y 
Li 
2 hiv 
or (CLM + Ca Wo) = EB [ (-v) sin k xdx = - 5 
Continuing, from equation (15) , 
5 2 2 
Vp 4 as dey 
Ox See, 2 oir ca 2 
Tt=0 GAK GAK 


But from equation (20), this means that, 








4 6 ; rts) . 
Jha on 
GAk GAk 
Nn n 


Since, as will be seen later in this work, w 


contradiction results, and 


C =O. 


2 amma er 


Lastly, apply equation (16), so that, 





Oo O 
3 Vp in 
xB i es in “in a hie Zz 

n 
G=o 
=. 


=e 


ata 


cannot equal Wo? 

2 

: Cc : Won 
3n Won 


a 


(20) 


(23) 








By combining equations (21) and (23), 


(Gak - - mw =) 

el aon en (2h) 
ln on ; “a Gn 0) 

La on liga! 

Co 

(GAk 2 - mw 
C = aad _ nn (25) 
sn Ty (w wo) 

on on in 


Thus, by applying equations (22), (24) and (25) to equation (12) and 
adding the quantity, vt, as indicated in Definition of Problem above, 


one finally obtains 


GAk 2 mw s+) 


2 ( 
4 n> min Sin w. t 
= = ol oe, ace te. 
ck) =ive2 ot nr rt, Cl en 


n=1,3,5... ft na 
) 
ce ee wt {| sin kx. (26) 


solution of the Frequency Equation. To "complete the solution efsrme 
differential equation of motion, the two series of frequencies must be 
evaluated, i.e. equation (11) must be satisfied. 

Thus solving equation (11) for we , one finds that 


EIm | + he EIm *|- JmEI , 4 
2 [nm + (oR + J) Sad - = (oR +d) "S * GA a 


; 21 -  —e 
> gm 
GA 


Clearly, there are two frequencies associated with each Kk, Po makblss, walerae 


is a two-fold infinity of frequencies, each pair of which is associated 


V2 





with an odd numbered mode of deflection. Intuitively, since the deflection 
was broken into components from shear and from bending, one would expect 
one of these frequency series to arise from the bending and the other 

from the shear. This point will be demonstrated more rigorously in a 

later section. 


Tt has been found more convenient to put equation (27) into the form 





ocak © 
W a 
“al 2 / 2 
G L IL G L iL LG 
mr n Tr n 
(27a ) 
or, letting 
2 
G L 1 
Wee n 
a 
D eGAk 





no ee oe) 
~ E 
Details of this derivation are given in Appendix B. 


Calculation of Shear and Bending Moment. Based on the equations derived 
in the previous sections, the values of shear and bending moment may now 
pescalculaced. 

Combining equations (22), (24), (25), and (17) the shear force at 


any section, at any time, is 


atc 





2 (cak © 
hi we in ) 
na=l5 5 Der. me er in 
(G in ° 
Ste 5 ' 2 w,, Sin w, t cos kx . (28 ) 
n eon pba 

















00 3 p 2 2 
mw) mW GAk - mw 
Ohf _ hv en + Jik w = ae ae eS Sin w~ t 
Ox } nt KO n en GAk mente am i 2 me 
Wendes oy tive en in 
mW, mu, > cak & mu) 
- kK. + J On ~ CAk 7“ 5 a 2 Sin Wn” cos k x 
en a 
Integrating this, 
= mW mw 3 GAk - mw = 
hiy on o n in 
M = =F el k W aos ee ae ie 
nttk kK nen GAk m (w a w,) {sin wt 
N= 3) San en as] en 
m mw GAk, © - mW, 
-|=— +4 K Wn > ane ae PS wt sin kx, (29) 
n n m (w, - Ww 
GA en in’ | 


where once again boundry conditions preclude a constant of integration. 


Evaluation of Equations. An indication of the validity of any equation 


occurs when the equation reduces to a simpler, classically known form 


St ye 





based on restrictions of the parameters. Accordingly, look at the values 
of y, Q, and M when L/r' goes te infinity. 


Thus, from equation (2Ta), 














W 2 eo 

Da “> ’ 
2 2 eas 

me BAK [rT 16 7 ic EI 

ln ? > nm . 
mL 

or 
a ye eee (30 a,b) 
on eiilia n m : 


Notice that Wn is the classical natural frequency for beams of long 


length to radius ratio. Putting equations (30 a,b) into (26) then 





& 
Lim _ y; ie Ale 
y (x,t) =ve t - noe |e Sin w, t sin kx 
Ly » © M=1 ,3 eps. in n n 
r! = 
or, ~ 
(x,t) =v t WBE fe sin x2 [Ee sin k (31) 
OS f n ne i en amide as Co 
iQil AGS Sp tarae 


But this is just the classical deflection equation for large =, [2] ! 
Now, look at the bending moment equation, equation (29). At first 


glance, it appears that the coefficient of the Sin w ae term blows up 


2 


-15- 


——_ 
_ - 
. | 
—- «- - 
- = - 





as — o , However, it can be shown that this term actually goes to 
zero. This is because the mass moment of inertia, J, times kK goes to 
zero. see Appendix C for proof of this convergence. 


Appendix C also indicates that 





C 
. mw 
ae M = av ue Sin W t Sin keke, 
7 n Ke din n 
Tare saa: “aia Nail 55. n 
16 


M a » Bim sin k* jELy . sinkx . (32) 
TT n n m n 


To further check this, it can be seen by inspection that using equation 


(31) to determine M from the formula, 


equation (32) results. 
The shear equation does not converge when = o , as is the case 


when classical theory is utilized. This fact shall not be derived here. 





Next, look at the equations for y, Q, and M when = goes to zero. 
Notice that this can occur in a beam of small =, for any mode, ns) Gmener 
any beam at a large n. 

Substituting equation (27b) in (26), (28), and (29), it can be shown 


that 


~16~= 





Gilt 





! - 2 
LI, [7m i iG 
y( xX ; 40, ) —_ V 16 a — —n =. = a _ ate | he lellenemene 





_ i 
esecee ya sin Ww, t - i - 
eo ee 2 25 aie 
E Z-,/z2 -— 
i 
ea L 
2 ke e 
Ze Z a a 
scouees > el sin Wt sin kx ) (33) 
7, = oes 
iE 
© Cc hia 
a 1 Z+i/Z Saar a 
a(x,t) = yah : 4 \——_=— 
n=1,3,5 
L 
eC 
ee Re! 2 he 
a(n - vo" - 2 Fer Ese ee ee 2 
coeee 2 ie sin w,t - i 
E 
i 
oO 
2 (2 + fe - ) 
eee - sin wt cos kX , (34) 
Phe 
EK 
& 
2 





n=1,3,5.. = eae E 
L 
ee a 
(2 + ve" - 2) 
+ e( wy sin 0% sin kx : (os) 
2 mas 


Be 


Certainly these equations look forbidding, but their derivation involves 
only algebra, even though long and complicated, and shall not be indicated 
heres 

Taking the limits of (33), (34), and (35) as L/r' goes to zero, one 


finds 


Lim hy, fa 
y(x,t) ac "2 GA 


© 

i 
la wat sin KX » (36) 
at eres 


L n= 
ri> 0 - 
© 
Lim hv /mGA He 
Q(x,t) = + = sinw, tcoskx , (37) 
n 
=> O n=1,3,5 7 
eo G 
lim M(x,t) = Avi yeh _ e sin wt 
th (2 ~ iE 
=> 0 Ma, 3,26 EK 
penne et Sines or. (38 ) 
@-) et] 
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The frequencies in the above equations are found by 


2 
lim 2 CGA, 


nay MEDEE DB) 


from which one sees that 


EA © GAk © 
» n w. © 2 (39 a,b) 
on? “mn ‘ In? mn ‘ 








One might ask why such complicated equations as (33), (34), and (35) 


were used to obtain the limits instead of merely substituting (39 a,b) 


“Whe 





1 





into (26), (28), and (29). Had this been done, the sin wt portion 


of the moment equation would have gone to zero, yielding a faulty 


equation. 
Equations (39 a,b) indicate that w,, is associated with the bending 
of the beam since it contains Young's modulus. Similarly, Wn is asso- 


ciated with the shearing action. This holds true for all modes in beams 
with small L/r' and for any beam in the higher modes. One could con- 
sider that the discrepancy in the lower modes is caused by the inter- 
action of shear and bending. There is not sufficient time for the 
modification to occur in the higher modes. 


With this interpretation, it appears that the response of a very 


short beam is almost entirely due to shear, which is completely acceptable 


intuitively. Only the bending moment shows any effects, then, of rotary 
inertia, since only it has the sin W,% term. 

A zero-length beam probably has no significance physically, but, 
when it is recalled that as n gets large, the nth term goes to the form 
of (36), (37), and (38) for all beams, they become more interesting. 


Notice equation (37) in particular. This equation can be transformed by 


fee) fee] 
hyo oa 1. —— uf 
a mGA A n ( sin wo 8 eos k x) B. hy y mGA n oe 
Mil seautoiees 2 Ww n=1,3,;50—0% 


(sin w, _t cos k x + cos w_ t sin k x) - eos bot sin kee 
ln n n n Es ln n 


Prope | 


i 
io] 
: bey 7 ial 4 /GA 
- sin Wnt cos 2) | ee mGA = [5 sin K (x $ a 
N=1,3,D-ce 


21o2 


5 sink, (x [SA 9) : (40) 





= =) 





This is just the equation for square waves traveling in from the ends. 


The magnitude of each component is v a The component traveling 
in from x = 0 is initially positive, and the one from x = L is initially 


Nefavaye. The waves travel at a velocity, 


velocity = (= = iE 


Further notice that, after each component has traveled the length of 
the beam once, a reflected shear wave of magnitude v /mGA results. 
Thus, the transient response of a shaft could be approximated by this 
Square wave, and the accuracy improves as =, gets small. 

Lastly, it is noticed that equations (36), (37), and (38) prove 
the convergence of the basic equations for y, Q, and M, since, as 
indicated above the terms of each of the equations represent the Limit 
as n gets large. Thus for the higher modes, Ve and Ms decrease as 
*, and om decreases as = 
Concept of an Equivalent Beam. An important characteristic of the 
E(uavlons Hor caetlecGion, shear, and bending moment can be realizedeaint 
they are non-dimensionalized. 


Accordingly, define non-dimensional values with primes. Then looking 


at equation (27b), one finds that 
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The non-dimensional time is defined such that 


wt ti = wt ; 
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t — = 
x = 1 (43) 


Continuing, from equations (33), (34), and (35), it is seen that 
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From these equations, it is clear that the non-dimensional values 
of frequency, deflection, shear, and bending moment are equal for any 
non-dimensional time or position, for any two beams if their values of 
Z are the same. 


Recalling that 


ee 


eal kop 


2 
iP 1 
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val 


it is thus necessary that, for the non-dimensional equations of any two 
beams to be equivalent, the ratios = and =, for one beam must equal 
these ratios for the other. 


Look at the ratio of the modified shear modulus to Young's modulus. 


As was indicated earlier, 
G=k'@G@ , 
where k' depends on Section shape. Evaluation of k' is complex and beyond 


the scope of this work; however, for a first approximation k' could be 


evaluated by the equation, 


2 
lohan 
Koo = 2E LC, (47) 
(° naa 
O 
where 
b = width of the beam material at the neutral axis, 
C = height of the extreme fiber of a section measured from the 
Reublral, oxXiLs, 
h = height of any point in a section measured from the neutral axis. [4] 
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See reference [5] for the exact evaluation of k'. Use of equation (47) 
yields results accurate to about four per cent for a solid circular cross 
section. 


Clearly there are many variations in shape that would yield a 
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constant value of 


Effect of Damping. Even though the equations derived herein are appli- 

cable only in the elastic region, hysteresis damping does occur. If 

the equations for y, Q, and M are used in the form of equations (26), 

(28), and (29), respectively, there is a possibility that many of the 

modes will become additive at some relatively long time. The result would 

be inordinately large magnitudes for the three quantities of interest. 
Applying damping in an analytical way would be difficult to say the 

ieasuesmiowever, 1ts major effect, i.e. that of reducing amplitude; @cemla 

be indicated approximately without loss of generality for the formulas. 
Accordingly Let the deflection energy that is lost each cycle be 

two per cent of the energy of that cycle [6] . From this, it can be 


shown that a damping factor for each term would be approximately 
Wee 
2 n 
Damping factor =e 0. 00161 (48 ) 
nee Appendix D for details of this derivation. 
cuggested Design Method. A method of calculating the response of any 


beam is implied by the concept of an equivalent beam. Using equations 


(42), (44), (45), and (46), modified by damping factors, one could 
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obtain values of maximum shear, deflection and bending moments and the 
time at which these occur. For a constant : » these maxima could be 


plotted versus =, - several plots for various = would result ina 
set of curves that would be completely general. After a few simple 
calculations to convert from the non-dimensional magnitudes to absolute 


magnitudes, one would have a good indication of the stresses experienced 


by the shaft of interest. 
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fit. RESULTS 


The equations for deflection, shear force and bending moment, 
modified by the appropriate damping factors, were programmed on the 
IBM 090 computer for a solid steel, circular beam. The results of 
that computation are indicated in Figures II, III, and IV. These plots 
are for unity excitation velocity; thus one need only multiply the 
actual velocity to obtain the magnitudes desired. 

Figure II shows the time history of y, Q and M for a beam of six 
inch radius and thirty-six inch length. Notice that, for a unit velocity 
input, the deflection curve magnitude could be considered in units of 
time. The origin of each shear and moment curve is the horizontal 
line representing the position of the end points on the corresponding 
deflection curve. Since the excitation is symmetrical only one-half of 
the beam is represented. 

Figure III is included to obtain a qualitative picture of the com- 
parison of the responses of two beams of different =, . the solidibine 
is for the beam of Figure II while the dashed line represents the response 
of a beam of six inch radius and sixty inch length. The units of this 
figure have been normalized by dividing time by length and y, Q, and M 
by the product of velocity and length. Since, for both beams, the cross 
section is the same, equations (42), (44), (45), and (46) indicate the 
validity of this normalization. Notice that the magnitudes have not been 
non-dimensionalized, although that too could have been done. Nevertheless 


the normalizing method employed provides a satisfactory base for comparison. 
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Figure IV shows the variation of the magnitude of the maximum 
deflection, shear, and bending with respect to =, . Here all units 
have been non-dimensionalized in accordance with equations (44), (45), 
and (46). The position of the maximum deflection and bending moment 
was, aS would be expected, the center of the beam. The position of the 


maximum shear varied from x' = 0 to about x' = .125. 
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Figure IV 
Non-Dimensional Maximum Deflection, Shear 


Force, and Bending Moment 
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IV. DISCUSSION OF RESULTS 


Upon examination of Figure II, one's attention is immediately drawn 
to the negative motion that precedes the positive deflection as the 
disturbance travels into the beam. The curves for shear and bending 
moment tend to support this, but no physical reasoning is evident. To 
show how this occurs, look at the solutions for ‘e and ui By putting 
equations (24) and (25) into (18) and (19), it is seen that, relative to 


the end points, 
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Equations (49) and (50) show that the coefficients of each term in 


the series is the same as the coefficient of the corresponding term for jy, 
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the net deflection, multiplied by a factor. For simplicity the nth term 


of equations (49) and (50) could be written, 
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where (1-) indicates that this factor approaches one. Now, by examining 

equation (26), the equation for y, it is seen that the coefficient of the 
sin w4° term is negative. Furthermore, - is always greater than one, 

e.g. for a solid circular shaft it is approximately 3.45. With these 


facts it can be seen that, for a circular shaft, 


ee wee S| Y, > (53) 


Vicia 4 _ 245 an . (54) 
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In other words, this shows that the deflection due to shear is greater 

than the net deflection and that the deflection due to bending is initially 
in the opposite direction. When it is recalled that the higher series of 
frequencies is associated with bending the the lower with shear, it fol- 
lows that the bending effects should travel into the beam faster than 


the shear effects. The ratio of these speeds will be on the order of 


Utilizing these concepts, Figure V shows, in qualitative terms, how 


the resultant deflection occurs. 


Figure V 
Qualitative Description of Deflection Components 
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fiecat be deaucen from this @iscuscsion, then, that the negative shear 
that precedes the expected positive values is due to the negative bending, 
and not the impact shear. This accounts for the fact that the negative 
portion keeps pace with the bending action instead of moving at the 
Slower speed of the shear wave. Conversely, the rising, poSitive bending 
moment 1s due to the expected positive shear. 

By equations (53) and (54), it can be seen that the energy contained 


in &@ beam is much greater than is indicated by the deflection alone. Even 


Be 





though the deflections due to shear and bending are of opposite sign, their 

respective energies are additive. Thus the actual energy in a beam is 

much greater than classical theory predicts. Lack of time precludes 

deeper investigation of this point to gain quantitative information. 
iuevaeenovrced CoaGiehe curves for shear are relatively Udo wane 

ragged. Even though the computations went to the forty-ninth mode, this 

is probably not enough since the series converges only as = , and as that 

Omi ye an bie higher modes. 

Upon examining Figure III, the lack of similarity in the two responses 
is noteworthy. This confirms that beams of two different length to radius 
of gyration ratios cannot be compared, even in a normalized form. Fur- 
ther, it is noted that, as predicted, the shear curves of the beam with 
tne longer =, look less like the ideal square wave than those of the 
ehormser beam. 

The accuracy of the curves in Figure IV is doubtful. Because of com- 
puter time limitations the search for the maxima may have been done with 
the time and position increments too large. For this reason also, the 
range of time examined may not have included the absolute maxima. However, 
the curves are included to show the general variation of the indicated 
quantities with =, ° 

The design method suggested in the Procedures section is dependent 
upon the availability of sufficient data from parent models. To obtain 
sufficient data to be generally useful, machine computation will probably 
be necessary. However, programming the equations derived in this report 
would not be too difficult. An example of how to utilize such derived 


curves is indicated in Appendix E. 
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A further refinement of the results of this work could be accomplished 
by incorporating them into elasticity theory to obtain the actual maxi- 
mum stresses experienced in a beam. A calculation of this sort could be 
included in the computer program previously mentioned and a single family 
of curves of maximum stress versus =, for a range oF 2 would result. 
Lastly, it is recognized that the results described herein should be 


confirmed by experiment before any actual engineering use of this is 


attempted. Such experiments are therefore recommended. 
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V. CONCLUSIONS 


It is felt that the equations derived speak for themselves and require 
no further qualification. The equations show that the deflection due to 
shear predominates in the total deflection and that the deflection due to 
bending is in the opposite direction. The energy contained in a beam is 
greater than is indicated by the net deflection because of this fact. 

The equations imply that the disturbing effects of bending travel through 
the beam faster than those of shear. 

It has been shown that beam response may be placed in a non-dimensional 
form whereby all beams with the same ratio of length to radius of gyration 
and the same ratio of shear modulus to modulus of elasticity, modified by 
a factor accounting for section shape, will have identical non-dimensional 
responses. 

This fact forms the basis of a simplified method to obtain design 
calculations. The method is, however, dependent upon the compilation of 


response data from a number of parent models. 
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VI. RECOMMENDATIONS 


It is recommended that the equations derived in this work be 
confirmed by experiment. 

The indication that significantly greater energies are contained 
in a beam than is implied by the net deflection should be investigated 
farce x . 

It is recommended that graphs for the non-dimensional form of the 
response of any beam be obtained. In doing so, care should be taken to 
assure that enough terms in the series be included to provide sufficient 
accuracy and to assure that enough of a span of time is examined, so that 
certainty of obtaining the maxima results. 

Lastly, it is recommended that the results obtained in this work 
be applied to the theory of elasticity so that prediction of the magnitude, 
time and position of maximum stress could be calculated. It is recommended 
that the results of this be placed into a family of curves, where each 
curve is for a constant modified shear modulus to Young's modulus ratio, 
and the coordinates are maximum stress and the ratio of length to radius 


of gyration. 
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VII. APPENDICES 
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Appendix A. Application of Boundary Conditions 


Equation (8) specifies that X(O) = X(L) = 0. However this could 


be specified more completely by, 


X, (0) =X, (Lb) =X, (0) =X, (L)=0. 


Now let C, = CLs + Ch, ; C, = Ca, + Coy? and so on. From this, one can look 


at the boundary conditions for the deflection due to bending separately. Thus, 


7h ee 2 A pe 
ron) C,, Sin kx + C,, cos kx + C3 sirh ‘kx + C), cosh.ikx , (a) 


subject to boundary conditions 








0 aaa (b ) 


H, (0) = Co, + Sy, = O » 


Therefore Con = Ch, =i: 


Then 


x, (L) = C,,, sin KL + C2 simh ‘kL =O , 
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ano 1b 3b 
x=L 
since sin hkL can never be zero, the above can only be true if Cop = 0 and 
he—onn 4 fh = 1, 2, or eooee 
Thus , 
: nit 
x(x) = C), sin > x. (c) 


But the constants for Xo the deflection due to shear, have not yet been 


satisfied. Thus, 


Xe (x) = C,, Sin kx + C,. cos kx + C25 sinhikx + C) , coshikx , (ad) 
and X. (j= Co, t+ Cy, = 0, 

X. (L) = oe sin kL + C,. cos kL + C3. Sith kL + C) . cosh kL =O. 
But from (ec) k = a one 

X. (L) = Cy, + C3. sinh nt + C) cosh.nn 

= C,. (1 - cosh nr) + C2. sinh nt = 0. 
But this is only true if C, =C,. = 0. Thus, X = (C__ + C,,) sin == 
2s 3s 1s Lb L 
= C sin = 


Since the excitation of the beam motion is symmetrical, the slope of the 


beam at x = 2 must be zero. This precludes any of the even numbered modes, and 


x = C. sin ES ae lion. + 3-6 (e) 
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Appendix B. More Convenient Form for Frequencies 


Starting with equation (27), multiply numerator and denominator by, 
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heculoing in, 
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For the case under study, a solid circular beam, = =— 


Thus, 
h 
(2 2EIk, 
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Now multiply numerator and denominator of the last term under the radical 


by at to obtain 
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Thus, by putting this into the equation for Ww » and factoring, 


wo = ————_4 = (b) 
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Appendix C. Proof of Convergance of Bending Moment Equation for Large 


Length to Radius Ratio. 


From equations (29) and (30 a,b) , it can be seen that, letting NT 


be the nth term in the infinite series of M, 
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Remembering that, 
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it can be seen that (a) goes to the form > as =. goes to infinity. 
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Clearly L' Hospital's rule applies here, but, in order to keep the 
application as simple as possible, recall that this calculation is only 
to show that this limit goes to zero. As the numerator is differentiated 
in accordance with L'Hospital's rule, it will always have zero as a limit 
as =, goes to infinity. Thus the case will be proved simply by showing 


that the limit of some derivative of the denominator becomes other than 


zero. 
To proceed, then, 
2 
2 of (L\ b 
Ge roe EA? r 
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2 a8) 2 i: ee 
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oince we are not interested in the actual magnitude of the limit, cancel 


out factors which are constant with respect to =. Thus we get 
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Therefore the sin 56 term goes to zero in the limit as =, goes 


LOwiLNiena ty. 
By similar reasoning, simply by multiplying J by 2 in the sin W, 4% 
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Appendix D. Derivation of Damping Factor 


A beam experiencing transverse motion could be interpreted as an 
extremely complex spring-mass system. Since the energy of any spring- 
mass system is proportional to the square of its deflection, one could 


write 
2 
where 


H, = Energy of the ith cycle, 


¥. = WeeleeLion of Ghe ith cycle. 


If the energy lost per cycle is two per cent of the energy of that 


Cycle. when 


E y = yo y 
ee aE itl 3 = = 0.02 . (a) 

= Y 

i “a i 


From elementary vibration theory the damping factor is of the form, 


-w Bt 
Damping factor =e ” , (b) 


where 
B = fraction of critical damping, 


and the frequency, modified by damping is 


2 
w, = wo Vi-B . (c) 


From (c) the time for one cycle, or the period,is 


7? 2 oe (d) 
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Cleary, theme crom (b) and (a) , 
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Combining (a) and (e) one obtains 





2 
y 
eee, -'nB a 
2 
V5 
Solving this for B one obtains 
B = 0.0@161 . (f) 


Notice that the effect of B on the natural frequency in (c) above 
is small. Therefore this effect shall be neglected. 

If each term in the equations for jy, Jig and Vy is multiplied by 
its appropriate damping factor and the equations for Q and M are derived 
again, it can be seen that, since B<<1 , the same equations result 
approximately except that each term is multiplied by the appropriate 


damping factor. 
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Appendix E. Example Problem 


Assume we have a solid steel circular shaft 18 inches in diameter. 
What would be the maximum length that could be allowed between supports 
if the beam is not to deflect more than .1 inch when subjected to a shock 
velocity of 35 inches per second? 


Let steel have the characteristics, 


oO 
- 7,33 x Sohail 1lb.-sec. 


ans 


g' = 11.2 x 10° psi. 


It can be seen that, by equation (47), 


KI 1} a = Bex Fa 5 
ae ey < aa 2 
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Oo 
Thus 


G= ? G' = 8.45 x 10° psi. 


By equation (44), 
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For a circular cross section, r' = x 





8.45 x 10° 


7.33 x 107! 






Tabulate y' and =, for several values of L, as indicated below: 


L (inches ) y' (35 =, (5) 
96 cial 21.3 
120 2.56 26r ii 
144 2645 32.0 
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Figure VI 
Determination of Maximum Length of Beam 


According to Example Problem 
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Plotting these values in FigureVI,.which is just a copy of Figure IV, 


it can be seen that the required L/r' is 30:3.) (hen 
L 
L=<=,xr' = 30.3x 4.5 = 136 inches. 


The maximum bending moment for a beam with this L/r' is obtained 


from FigureVI and equation 46, i.e., 


M' (vL VmGA ) = M' (vLAV 9 ) 
35 x 136 xm7x 9° V7.33 x 107" x 8.45 x 10° 


419,000 in.-lbs. 


= 
I 


Using the formula for maximum bending stress, 


Db. ee 2 
we find 
S. = a T23mpsi.. 
Ee) 


Similarly, from equation (45) and Figure VI, 
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OF=—0. (v /mGA ) =. (5 Keene TX 9° 7.33 x 10 x 8.45 x 10 = 2010 ib 


Now, the maximum shear stress is 





=—ElOese DSi, . 


Clearly the deflection governs in this case, but if the allowable 
deflection had been greater, the acceptable length would have increased. 
This would cause a rise in bending stress. A point would thus be reached 


where the stresses govern the length. 
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